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IS  AUJTHAC  T 


The  first-order  statistics  of  terrain  clutter  are  investigated  theoretically  and  experimentally. 
The  return  from  a  point-scatterer  (a  scintillating  signal)  is  used  as  a  basis  for  generating  the 
statistical  properties  of  terrain  clutter  as  observed  by  side-looking  airborne  radars.  Accord¬ 
ing  to  the  model  the  limiting  distributions  of  terrain  clutter  are  the  Gaussian  distribution  and 
the  distribution  for  a  scintillating  signal  (specular  clutter).  An  array  of  point-scatterers  within 
the  radar  beam  produces  clutter  which  has  long  tails  similar  to  those  of  a  lognormal  distri¬ 
bution. 

The  convergence  of  the  first-order  statistics  of  a  scintillating  signal  to  the  statistics  of  a 
sinusoidal  signal  is  demonstrated.  This  property  allows  the  use  of  many  statistical  results  avail¬ 
able  in  statistical  communication  for  sine  waves  in  Gaussian  noise  to  interpret  and  predict 
in  many  cases  the  statistical  properties  of  terrain  clutter. 

Also,  the  statistics  of  terrain  clutter  are  inferred  from  data  taken  with  the  Naval  Research 
Laboratory’s  Four  Frequency  Radar  System  by  the  Kolmogorov— Smirnov  test  of  the  cumula¬ 
tive  distribution  and  by  the  computation  of  the  first  five  central  moments  of  the  normalized 
radar  cross  section  of  terrain  (in  decibels).  It  is  found  that  in  general,  terrain  clutter  statistics 
is  not  exponential  (with  Rayleigh  envelope)  nor  lognormal  distributed,  which  supports 
some  of  the  predictions  of  the  analytical  results. 
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POINT-SCATTERER  FORMULATION  OF  TERRAIN  CLUTTER  STATISTICS 


INTRODUCTION 

Through  the  years  terrain  clutter  has  remained  as  the  main  limitation  of  airborne 
radars  in  the  detection  of  targets  when  operating  over  land.  To  describe  terrain  clutter  in 
a  more  realistic  manner  and  help  radar  designers  and  analysts  to  make  provisions  for  the 
large  probability  of  false  alarms  obtained  in  this  environment,  many  attempts  have  been 
made  to  understand  the  basic  phenomena  responsible  for  this  clutter  and  then  be  in  a  posi¬ 
tion  to  predict  the  effects  of  terrain  clutter  on  the  radar  systems. 

Presently,  it  is  pretty  well  accepted  that  terrain  clutter  is  the  result  of  two  basic 
mechanisms  (1,2):  strong  point-scatterers  (or  speculars)  and  a  Gaussian-distributed  clutter 
of  Rayleigh-distributed  envelope,  which  is  produced  by  many  equal-amplitude 
scatterers. 

In  previous  statistical  investigations,  a  priori  statistics  (e.g.,  lognormal,  Rayleigh,  Chi- 
Square  and  Rician)  known  to  occur  in  terrain  clutter  have  been  used  from  the  outset  to 
describe  the  process.  Here,  a  much  more  fundamental  point  of  view  has  been  taken.  The 
complete  process  is  developed  in  terms  of  the  statistical  properties  of  the  scintillating  re¬ 
turn  from  the  elementary  point-scatterer,  and  the  more  complicated  distributions  can  be 
obtained  by  a  superposition  of  many  point-scatterers  within  the  radar  beam.  Although 
the  basic  idea  is  quite  general,  in  this  investigation  it  will  be  applied  for  the  particular 
case  of  a  side-looking  radar. 

With  the  point-scatterer  as  the  building  block  of  the  clutter  process,  a  great  deal  of 
physical  insight  can  be  obtained  and  terrain  clutter  data  may  be  interpreted  from  a  more 
fundamental  point  of  view.  For  example,  point-scatterers  in  many  instances  can  be  identi¬ 
fied  in  a  topographic  map  or  in  a  more  precise  manner  by  means  of  synthetic  aperture 
maps  obtained  at  the  same  radar  frequency.  That  is,  they  are  in  most  cases  physical  en¬ 
tities  which  appear  In  the  form  of  various  natural  terrain  features  or  man-made  objects. 

In  the  analytical  investigation  it  is  proved  that  the  first-order  statistics  of  a  scintillating 
signal  (i.e.,  a  linear  frequency-modulated  signal  produced  by  a  point-scatterer  as  the  air¬ 
borne  radar  flies  by)  converge  asymptotically  to  the  first-order  statistics  of  a  sinusoidal 
signal,  or  more  precisely  they  converge  in  density  and  distribution  to  the  statistics  of  a 
sinusoid  for  large  frequency  modulation  rates  and  large  observation  times. 

Analytical  expressions  are  derived  for  the  first-order  statistics  of  a  scintillating  signal 
using  the  assumption  that  its  frequency  spectrum  is  constant  within  the  time  of  the  fre¬ 
quency  sweep.  In  the  more  realistic  case  in  which  several  point-scatterers  appear  within 
the  resolution  cell,  it  is  not  possible  to  obtain  explicit  analytical  expressions  for  the  first- 
order  statistics  of  the  return.  In  those  cases,  the  statistics  may  be  obtained  by  convenient 
computer  simulation  or  by  invoking  the  asymptotic  property  of  the  statistics  of  a  scintil¬ 
lating  signal.  In  many  cases  it  is  possible  to  obtain  v  n  •  v.  dr  l  of  insight  by  proper  inter¬ 
pretation  of  available  statistical  results  for  sine  waves  anu  Gaussian  noise.  Many  of  the 
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statistical  properties  for  the  envelope  of  terrain  clutter  may  be  obtained  in  a  similar 
manner. 

The  convergence  of  the  statistics  of  the  sum  of  n  sine  waves,  of  irrational  periods  of 
the  same  amplitude,  to  the  Gaussian  distribution  is  of  order  1/n  as  n  -*  °°.  However,  the 
convergence  of  the  statistics  of  n  scintillating  signals  of  equal  amplitude  should  be  of  order 
1  /y/n  or  slower,  since  the  mean  value  of  scintillating  signals  is  nonzero  and  the  component 
signals  are  not  independent.  The  slower  convergence  of  the  statistics  of  the  sum  of  n 
scintillating  signals  may  be  directly  responsible  for  the  lognormal  statistical  properties  of 
clutter  which  sometimes  are  observed  in  terrain  clutter. 

Finally,  the  statistics  of  terrain  clutter  obtained  with  the  Four  Frequency  Radar  (4 
FR)  system  are  obtained.  The  empirical  investigation  shows  that  in  general  terrain  clutter 
is  neither  Gaussian  (with  an  exponentially  distributed  radar  cross  section)  nor  lognormal. 
These  conclusions  have  been  arrived  at  by  means  of  a  Kolmogorov-Smimov  test  of  the 
cumulative  distribution  of  the  normalized  radar  cross  section  (NRC)  of  terrain  and  by  an 
investigation  of  the  first  five  central  moments  of  the  NRC.  As  a  general  rule  the  type  of 
terrain  is  reflected  most  markedly  in  the  mean  value  of  the  NRC. 


STATISTICAL  PROPERTIES  OF  SCINTILLATING  SIGNALS 
Scattering  Considerations 

Consider  an  airborne,  side-looking,  coherent  pulsed  radar  system  illuminating  the  ter¬ 
rain  at  small  depression  angles;  see  Fig.  1.  In  such  a  scattering  configuration,  the  range 
resolution  is  determined  by  the  radar  pulse  width  and  the  azimuthal  resolution  is  a  func¬ 
tion  of  the  horizontal  dimensions  of  the  antenna  (3). 

The  radar  transmits  a  characteristics  spectrum  of  electromagnetic  plane  waves  which 
are  scattered  by  the  terrain  in  all  directions.  The  backscattered  waves  are  collected  by  tne 
radar  receiver  and  they  constitute  the  terrain  clutter,  whose  statistical  properties  we  are 
interested  in  determining. 

We  will  be  dealing  with  the  simplest  form  as  possible  of  the  backscattered  signal  in 
order  to  include  only  the  necessary  features  of  the  return  containing  information  on  the 
physical  process  responsible  for  terrain  clutter.  For  our  purpose,  the  return  from  a  point- 
scatterer  within  the  resolution  cell  will  be  of  the  complex  fc'm 

Y(t)  “  for  0  <  t  <  T,  (1) 

keeping  in  mind  that  wt  are  interested  in  the  real  part  of  y(t).  The  continuous  wave 
(CW)  version  of  the  back  scattered  signal  about  the  angular  frequency  u>o  will  be  used  for 
simplicity.  A(t)  contains  both  the  amplitude  weight  of  the  two-way  antenna  radiation  pat¬ 
tern  and  the  intrinsic  amplitude  of  the  scatterer,  R(t)  «  \/(Rq  +  is  the  slant  range 
to  tne  scatterer,  <t>(t)  °  2 (}R(t)  is  the  phase,  v  is  the  speed  of  the  aircraft,  0  =  2ir/X  i"  the 
propagation  constant  of  a  plane  wave  transmitted  by  the  radar,  and  T  is  the  observation 
time  dependent  on  the  antenna  beamwidth.  In  many  practical  cases  vt  «  R o,  R o  being 
the  slant  range  from  the  radar  to  the  point-scatterer  at  broadside.  The  return  signal  for 
these  cases  can  be  approximated  by 
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Fig.  1— Scattering  geometry  for  a  aide-looking  radar 


Y(t)  =  €M0+K‘2)e->u0t  ,  (2) 

«0 

where  <t> o  ~  20f?o  K  -  (Sv^/Rq,  which  is  the  scintillating  rate.  For  t  <  0,  Eqs.  (1)  and 
(2)  will  apply  if  we  let  <p(t)  =  -  <p(~t),  since  the  rate  of  change  of  slantrange  in  this  case 
is  negative.  For  our  purpose  we  will  deal  with  the  return  for  t  >  0  only,  since  by  itself  it 
contains  all  the  intrinsic  properties  of  the  return. 

Thus,  terrain  clutter  may  be  represented  as  the  superposition  of  a  great  number  of 
scintillating  signals  at  any  given  time 

n 

Y(t)  -  for  0  <  t  <  T,  (3) 

*- J 

and  for  simplicity  ail  amplitude  effects  have  been  lumped  into  the  a*  coefficients  which 
are  taken  to  be  constants. 

Thus,  in  fact,  tacitly  we  are  assuming  that  the  illumination  pattern  of  the  radar  an¬ 
tenna  is  independent  of  angle  of  uniform  weight.  The  scintillating  rate  K  and  the  phases 
4>i,  are  also  assumed  to  be  constant.  This  is  really  a  simplification  because  in  practice  the 
amplitude,  the  phases,  and  even  the  number  of  scatterers  producing  the  return  may  change 
within  tire  observation  time  T. 
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One  Scintillating  Signal 

In  what  follows  we  will  deal  with  a  coherently  detected  version  of  the  scintillating 
signal  of  Eq.  (1).  The  fluctuations  in  this  case  appear  about  zero  frequency. 

Y(t)  =  A  cos  (Iff2  +<»  for  0  <  t  <  T,  (4) 

which  is  a  linear  frequency  modulated  signal.  The  frequency  spectrum  of  Eq.  (4)  is  a  com¬ 
plicated  expression  involving  Fresnel  integrals  (4). 

In  Fig.  2  the  frequency  spectrum  of  a  scintillating  signal,  as  obtained  by  Klauder  and 
others,  is  shown  as  a  function  of  the  parameter  KT 2.  As  KT2  -*■  <»,  the  spectrum  of  a 
scintillating  signal  tends  to  a  constant.  In  practice  the  side-looking  radars  operating  at 
small  depression  angles  KT 2  is  large.  In  this  analysis  we  will  asusme  that  the  frequency 
spectrum  of  the  scintillating  signal  from  a  point-scatterer  is  constant  during  the  observa¬ 
tion  time.  Thus,  from  elementary  probability  theory  the  characteristic  function  of  Eq.  (4) 
can  be  obtained  from 


where 


Fa)  = 


X  -  Kt 


p(X)  = 


1 

2*o 


1 

2 KT  ' 


Using  the  identity 


eiZcod  .  ^  imemJm  ( Z )  cos  md 
mm  0 


(5) 


(6) 


Fifc  2—  Frequency  *  pc  clrum  of  alineai  frequency 
modulated  tifnmi 
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where  eo  =  1  and  em  =  2  for  m  >  1,  in  Eq.  (5).  By  reversing  the  order  of  integration  and 
summation,  the  following  expression  can  be  obtained  for  the  characteristic  function: 


F(£)  =  2 ^imbmJm(lA), 

m«  0 


(7) 


with 


,  r“°  w 

Jm  2X0 


bm  =  f'  cos  ♦  <t>)dX. 


By  making  a  change  in  variables  it  is  simple  to  show  that 

/  7f  \  1/2  1 

bm  =  y  cos  m0C(rm)  -  sin  m<j>S(Tm) 


for  m  >  1 


and 


where 


*0  =  1  , 


with  C(rm )  and  S(rm )  being  the  Fresnel  integrals  defined  by 


and 


C(rm) 


S(rm) 


■I  “'(H1"' 


The  probability  density  function  (pdf)  is  by  definition  the  Fourier  transform  of  the 
characteristic  function.  After  the  manipulations  we  find  that 


p(Y) 


-  ir*1  (A2  -  F2r1/2  y  (-)"* bm  cos(m  sin’1  ~  ♦  ~  j. 

,  0  m“0  , 


m  <  a 
m  >  a 


w 


As  usual  from  the  pdf.  the  distribution  function  can  be  obtained  with  no  difficulty. 
It  on  be  shown  to  be 
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«■{»>*■> -(i-M)«  ^ 


m 


.  I  .  .1  V  mn\ 
mjmsin  *  Y)‘  (9) 


The  mean  value  of  the  scintillating  signal  can  be  shown  to  be 


,  ,  A  f  it  V2 

Mean{y}  =  — J  (costfC^)  -  sin05(ri)]  . 


(10) 


The  term  in  parentheses  outside  the  infinite  sum  in  Eq.  (9)  is  the  expression  for  a 
sinusoidal  signal  of  amplitude  A.  Since,  in  deriving  Eqs.  (7)  through  (10)  we  have  assumed 
that  KT 2  is  sufficiently  large  so  that  the  spectrum  of  the  scintillating  signal  is  constant  in 
the  interval  (0,  2 KT),  we  may  repla  e  the  Fresnel  integrals  by  their  asymptotic  value  0.5. 
Thus,  the  asymptotic  form  of  Eqs.  (8)  through  (10)  are  respectively. 


,  <4  2  —  v2\-l/2 

p(T)  ~  -(A2-T2)-1'2  ♦  - -j- 

»  2Ty/nK 


Z]bScos(f  +  m*) 

1 


(id 


v  [  .  »  y  mir\ 

X  cos^nsm-.-  ♦  TJ 


-»{v>v)~(i  - 


1  V-!  (-)"*1 


2 TsfiKL-i  m* 2 


cos^  * 


(12) 


sin^n 


v  .  i  ,  .J  mir\ 

X  sin  [  m  sin  1  —  ♦ 

A  j 


and 


(13) 


In  Fig.  3,  the  convergence  of  the  distribution  function  of  a  scintillating  signal  toward 
the  distribution  of  a  sinusoidal  signal  for  increasing  values  of  T  is  demonstrated.  The  cu¬ 
mulate  distribution  for  the  scintillating  signals  in  the  figure  were  obtained  by  computer 
simulation  by  sampling  Kq.  (4).  The  result  is  in  agreement  mth  the  prediction  of  Kq.  (12). 
For  example,  for  V*  *  0  and  0  *  0.  Kq.  (12)  assumes  the  form 


P{Y  >  9}  ~ 


1  \  ’  <-£•*»__ 

2 (2m-'l)3«  * 


04) 


and  the  numerical  value  of  the  infinite  sum  is  0.865  to  three  decimals. 


Y 

Kit.  3-TKc  distribotias  of  o tut  -ictaUlUtiat  us*il 

Thus,  for  large  7N/K  factors,  the  first-order  statistics  of  a  scintillating  signal  are  very 
wtiniir  to  those  for  a  sine  «w  atJid  as  T\fK  ■*  ®*  they  become  the  same.  However.  some 
important  differences  exist,  for  example  while  the  mean  value  for  a  sine  wave  is  aero,  the 
mean  value  for  a  scintillating  signal  is  nonzero,  given  by  Eq.  (17).  This  difference  in  the 
first-order  statistics  of  these  two  signs  it  yields  a  different  rate  of  convergence  oi  the  first 
order  statistics  for  the  sum  of  scintillating  signals  as  comjvared  to  that  for  sine  waves  (as 
live  number  of  live  components  increases),  toward  the  Gaussian  distribution.  This  point  is 
investigated  later  on. 

Central  Limit  Theorem 

As  an  illustrative  example  consider  the  convergence  of  the  distribution  of  the  sum  of 
n  sine  waves  with  irrabonal  periods  (irrational  periods  are  required  so  that  the  sine  waves 
are  mutually  indef  lenders  t )  and  amplitude*  c*  to  the  Gaussian  distrsbufron.  Via  the  charac¬ 
teristic  function  it  is  quite  easy  to  derive  the  pdf  which  u  of  ths  form  (S) 
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P(F) 


J: _  -y2/2a2J 


\/  2tw* 


L°> 

k=i 


64o4 


(f4  ~6f2  +3)  +  . . ; 


and  the  distribution  function  which  is  given  by 


Z>*4 


where 


(15) 


p{t>f}  ~  1 11 -*<*>]  +^^2,2jLe^-{3t~t3)  +  •■•’ 


(16) 


2a2  =J^ak\ 

k=l 


Y 

t  =  — 
a 


and 


<I>(f) 


■M 


eX2/2dX. 


The  convergence  of  the  density  function  for  the  sum  of  sine  waves  of  the  same  ampli¬ 
tude  toward  the  Gaussian  distribution  is  illustrated  in  Fig.  4. 


Fig.  4— The  density  function  of  the  sum 
of  sine  waves  of  e4ual  amplitude  and 
irrational  periods 


The  general  convergence  problem  of  the  distribution  of  the  sum  of  n  independent 
random  variables  to  the  Gaussian  distribution  was  investigated  by  Gnedenko  and  Kolmogorov 
f6).  For  the  more  general  case,  in  which  each  of  the  components  is  distributed  according 
to  a  different  density  function,  the  convergence  of  the  probability  density  function  of  the 
sum  to  the  Gaussian  distribution  is  in  the  following  manner: 
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vi ->'<••■*  l  ■**■-'>-< 


The  net  conclusion  obtained  from  examining  Eqs.  (15)  through  (18)  is  that  the  dis¬ 
tribution  of  the  sum  of  independent  random  variables  with  zero  mean  converges  to  the 
Gaussian  distribution  as  1/n  with  the  number  of  components  and  when  the  mean  value  of 
the  component  random  variables  is  nonzero  the  convergence  of  the  distribution  of  the  sum 
to  the  Gaussian  distribution  is  slower,  of  order  l/> fn~. 


I 


s 

j? 

5 

I 


Bernstein  (7)  investigated  the  convergence  of  the  distribution  of  the  sum  of  depend¬ 
ent  random  variables.  In  this  case  the  Central  Limit  Theorem  arguments  apply  to  the 
number  of  independent  subsets  contributing  to  the  sum,  which  obviously  is  smaller  than 
the  total  number  of  the  component  random  variables. 


Many  Scintillating  Signals 

In  general,  terrain  clutter  is  the  result  of  many  point-scatterers  within  the  antenna 
beam.  Analytical  expressions  for  the  first-order  statistics  for  the  sum  of  n  scintillating 
signals  of  the  form 


Y(t)  =  2^«fe«>s[A:(f-ffe)2  +  <pk] 

fe=i 


for  0  <  t  <  T 


(19) 


cannot  be  derived.  However,  several  important  properties  for  the  sum  of  many  scintil¬ 
lating  signals  can  readily  be  obtained.  For  example  using  Eq.  (13),  the  mean  value  of  Eq. 
(19)  is  immediately  obtained 


n 

Mean{y}  ~  ^  Cfe  c°s(j  +  <pk^  •  (20) 

ft=l 

Similarly  an  explicit  expression  for  the  variance  of  Eq.  (19)  can  be  obtained,  this  is 


where  5  is  the  Kronecker  delta  function  and  fai  -  (K/2)(tk  ~  t\  +  r)2  +  (4>k  +  0/)- 

It  is  not  difficult  to  arrive  at  the  result  that  the  worst  clutter  case  (greatest  proba¬ 
bility  of  false  alarm  for  a  given  threshold  level)  will  occur  when  all  tne  point-scatterers  are 
lumped  together  and  add  constructively.  For  this  case  the  variance  takes  the  form 
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and  the  first-order  statistics  are  those  for  a  single  scintillating  signal  of  amplitude  equal  to 
the  sum  of  the  component  amplitudes. 

Another  interesting  case  is  that  in  which  the  components’  scintillating  signals  can  be 
considered  mutually  independent.  In  this  case  the  variance  is  of  the  form 


n 


in  which  the  first  sum  is  the  variance  of  n  sine  waves  of  irrational  periods  and  of  amplitude 
identical  to  those  of  the  scintillating,  signals. 

Although  the  first-order  statistics  of  the  sum  of  several  scintillating  signals  cannot  be 
obtained  explicitly,  they  can  be  derived  very  conveniently  by  computer  simulation.  In 
Fig.  5  the  computer  simulation  has  been  performed  for  the  case  of  two-scintillating  sig¬ 
nals.  The  amplitude  of  one  component  is  kept  constant,  while  the  amplitude  of  the  other 
component  is  varied  from  0  to  1,  in  steps  of  0.2. 

Since  the  mean  value  of  scintillating  signals  is  nonzero,  the  distribution  of  the  sum  of 
many  scintillating  signals  should  converge  to  the  Gaussian  distribution  according  to  Eq.  (18) 
and  the  number  of  components  is  the  number  of  independent  subsets.  Here,  we  will  not 


,‘Ml  LITUDE 

Fig.  5— The  distribution  i unction  for  the  gum  of  two-scintillating  signals 
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attempt  to  define  what  we  mean  by  independent  subsets  in  terms  of  the  scintillating  sig¬ 
nals.  Our  aim  is  only  to  emphasize  the  slower  rate  of  convergence  of  the  distribution  of 
the  sum  of  scintillating  signals  toward  the  Gaussian  distribution.  This  is  illustrated  in  Fig. 

6  in  which  the  distribution  of  a  number  of  scintillating  signals,  for  point-scatterers  of  equal 
amplitude  in  an  array  configuration,  was  obtained  by  computer  simulation  for  different 
numbers  of  scintillating  signals.  The  limiting  distributions  are  the  distribution  for  one 
scintillating  signal  and  the  Gaussian  distribution.  Surprisingly  enough  the  distribution  for 
the  sum  of  a  large  number  of  scintillating  signals  has  long  tails,  similar  to  the  lognormal 
distribution.  Thus,  seemingly  we  have  discovered  a  new  mechanism  that  yields  lognormal 
clutter. 

Based  on  the  asymptotic  properties  of  the  first-order  statistics  of  the  sum  of  scintil¬ 
lating  signals,  the  distribution  of  the  return  from  a  point-scatterer  in  a  homogeneous  Gaus¬ 
sian  clutter  may  be  approximated  by  the  well-known  statistical  results  for  a  sine  wave  in 
Gaussian  noise  (8-10).  In  Fig.  7,  Rice  results  (9)  have  been  used  to  represent  the  distribu¬ 
tion  of  a  point-scatterer  in  various  amounts  of  Gaussian  clutter. 


I-PROBABILITY  DISTRIBUTION 
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Fig.  7— The  distribution  of  a  point  scaiterer  in  Gaussian  clutter 


In  terrain  clutter  measurements  taken  with  the  4FR  system  over  Phoenix,  Arizona, 
on  Oct.  7,  1966,  the  lognormal  and  specular  nature  of  clutter  was  quite  evident.  In  Fig. 

8  a  couple  of  these  interesting  distributions  of  the  complete  signal  are  displayed.  (The  out¬ 
put  of  the  amplitude  and  phase  channels  of  the  4FR  system  have  been  combined  to  form 
E  cos  <p\  E  is  the  envelope  of  the  clutter  and  4>  is  the  instantaneous  phase.)  The  distribu¬ 
tion  for  P-Band  exhibits  the  long  tail  typical  of  lognormal  clutter,  quite  similar  to  those 
given  in  Fig.  6  and  obtained  by  computer  simulation.  The  distribution  for  C’-Band  may  lx? 
a  composite,  produced  by  a  couple  of  strong  point-scatterers  in  a  background  of  weaker 
Gaussian  clutter. 


STATISTICS  OF  THE  ENVELOPE 


In  many  applications  one  is  interested  mainly  in  tire  statistics  of  the  envelope  E(t)  of 
terrain-in  particular  in  incoherent  radars.  The  envelope  of  Eq.  (3)  is  given  by 
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Fig.  8— Cumulative  distributions  obtained  with  the  4FR  system  ovor  Phoenix, 
Arizona,  in  1966  which  support  the  specular  nature  of  terrain  clutter 
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Interestingly  enough  the  mean  value  of  E2  is  given  by 


{£2}  =  ^P, 


for  <pk  *  <t>,, 


which  also  is  related  to  the  mean  power  of  the  incoherent  clutter. 

Thus  since  Eq  (26)  also  represents  the  mean  value  for  the  superposition  of  sine  waves 
of  amplitude  similar  to  the  scintillating  signals,  one  would  suspect  that  the  distribution  of 
the  envelope  of  terrain  clutter  may  converge  faster  to  the  Rayleigh  distribution  than  the 
distribution  of  the  complete  signal  (amplitude  ~”d  phase)  will  converge  to  the  Gaussian 
distribution.  Recalling  the  results  for  a  sine  wa.t-  Gaussian  noise,  a  point-scatterer  in 
Gaussian  clutter  should  have  an  envelope  which  is  Rician  distributed. 

The  convergence  of  the  first-order  statistics  of  the  envelope  of  the  sum  of  many  ran¬ 
dom  processes  toward  the  Rayleigh  distribution  is  reviewed  in  detail  by  Beckmann  (11). 

For  example  the  envelope  of  the  sum  of  n  sinusoidal  signals  of  irrational  periods  and 
amplitude  a *  converge  toward  the  Rayleigh  distribution  according  to  the  expression 


p(E)  ~  ^-f2/2°2{l  -  (8-8£2  + *4)  + 

0+  64  o* 


1 


for  the  density  and 


r\^>K 


for  the  distribution  function,  where 


t  * 
C  9  — 


iR*-2R2)  + 


Z>*3 


The  more  general  question  on  Uie  rate  of  convergence  of  the  envelope  of  the  sum  of 
independent  random  processes,  each  distributed  according  to  a  different  distribution,  can 
be  shown  to  converge  to  the  Rayleigh  distribution  according  to  the  expression 


fXE)  ~  1  +  i  (2-4*2, *4)  «. 
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for  the  density  with  £  =  E/M j^2  and 


f  E  ^  J  _  -*2  r  ,  1  (Af 4  -  2M<j ) 


V»i  >  1 


e  "  1  + 


for  the  distribution  function,  where 


M, 


4  Mi 


n 


fc=l 


M'2  = 

k-l 


(R*-2R2)  + 


(30) 


Ma 


!>«• 


and  m2fc  apH  /n4fc  are  the  second  and  fourth  moments  of  the  component  random  processes. 
Beckmann’s  expression  (38)  in  Ref.  11  is  incorrect;  a  quadratic  term  in  the  series  for  the 
distribution  is  missing.  Compare  with  our  (30). 


INFERENCE  OF  THE  STATISTICS  OF  TERRAIN  CLUTTER 

In  this  section,  the  statistics  of  terrain  clutter  us  observed  with  the  4FR  system  is  in¬ 
vestigated.  The  4FR  system  transmits  at  428  MHz  (UHF  or  P-Band),  1228  MHz  (L-Band), 
4455  MHz  (C-Band)  and  at  8910  MHz  (X-Band).  An  excellent  description  of  the  4FR  sys¬ 
tem  can  be  found  in  Ref.  1.  Since  the  amplitude  channel  is  completely  calibrated,  these 
data  will  be  used  in  the  inference  of  the  statistics  of  terrain  clutter.  The  outputs  of  the 
amplitude  channel  are  samples  of  radar  cross  section  of  terrain  (in  units  of  decibels)  as  a 
function  of  time  and  can  be  normalized  by  the  area  on  the  terrain  illuminated  by  the 
radar. 

The  data  used  in  this  investigation  were  collected  by  the  radar,  in  skU  looking  opera¬ 
tion,  on  Oct.  6,  1966.  over  Arizona  and  contain  return  from  deserts,  mountains,  farm¬ 
land,  and  local  communities.  Additional  information  on  tiie  measurements  can  be  found 
in  Ref.  12. 

Tile  cumulative  distribution  of  the  NUC  of  terrain  is  checked  with  the  Kotmogorov- 
Smimov  test  (13)  against  the  exponential  and  the  lognormal  distributions.  Finally,  the 
first  five  central  moments  of  the  samples  (in  units  of  decibels)  of  the  NRC  of  terrain  are 
obtained.  Clearly,  if  all  the  moments  were  known,  U»e  exact  distribution  of  terrain  clutter 
could  be  reconstructed. 

It  was  decided,  for  obvious  practical  reasons  of  terrain  identification,  that  the  total 
number  of  samples  used  in  each  distribution  and  moment  computation  should  not  exceed 
the  time  in  which  the  airplane  takes  to  travel  the  distance  equal  to  the  azimuth  antenna 
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beam  width  illumination  of  the  terrain.  Thus  the  total  number  of  samples  used  in  each  dis-  j 

tribution,  or  moments  computation,  is  a  function  of  airplane  altitude,  speed,  antenna  de-  j 

pression  angle,  and  antenna  azimuth  beamwidth.  The  total  number  of  samples  used  in  j 

each  computation  ranged  from  about  900  to  11,000.  j 


The  Distribution  Function 

The  Kolmogorov-Smimov  test  is  a  distribution-free  type  of  statistical  test  used  to  test 
the  hypothesis  if  the  cumulative  distribution  obtained  from  a  set  of  independent  samples 
taken  from  a  population  belongs  to  a  theorized  distribution  for  the  population.  The  num¬ 
ber  of  independent  samples  for  the  terrain  clutter  data  was  determined  by  a  standard  run 
test  (14)  on  groups  of  1024  consecutive  samples.  The  cumulative  distribution  of  terrain 
clutter  (between  the  0  and  the  100  percentile)  are  then  tested  against  the  exponential  and 
the  lognormal  distributions.  In  the  comparison  with  the  exponential,  the  mean  value  of 
the  theorized  distribution  was  adjusted  for  the  smallest  maximum  deviation  between  the 
distributions.  In  the  comparison  with  the  lognormal  distribution,  the  variance  and  the 
median  of  the  theorized  distribution  were  taken  to  be  that  of  the  sample. 

The  Kolmogorov-Smimov  test,  at  the  99%  level  of  significance  (i.e.,  the  cumulative 
distribution  is  accepted  as  belonging  to  the  theoretical  distribution  if  the  maximum  devia¬ 
tion  is  loss  than  or  equal  to  the  critical  value  as  given  by  Ref.  13)  shows  that  in  general 
the  NRC  of  terrain  is  not  exponentially  distributed  (with  Rayleigh  envelope)  nor  lognormal 
distributed.  However,  specific  sections  of  data  may  be  found  to  be  exponentially  distrib¬ 
uted  and  on  a  few  occasions,  for  small  independent  samples  sizes,  the  data  may  be  found 
to  be  lognormal. 

For  example  in  Figs.  9  and  10,  the  data  is  very  likely  to  be  exponentially  distributed, 
while  in  Figs.  11  through  15  the  data  are  neither  exponential  or  lognormal.  As  a  matter 
of  fact,  cumulative  distributions  which  have  long  tails,  typical  of  lognormal  distributions, 
fail  the  Kolmogorov-Smimov  test  in  comparison  with  a  lognormal  distribution  between  the 
0  and  the  100  percentiles.  Thus,  strictly  speaking  terrain  clutter  may  never  be  lognormal 
distributed  and  only  the  tails  of  the  distribution  may  be  approximate  by  this  distribution. 


The  Central  Moments 

For  purposes  of  identifying  the  statistics  of  terrain  clutter  in  a  more  definite  manner, 
the  first  five  central  moments  of  the  NRC  of  terrain  clutter  (in  units  of  decibels)  will  be 
derived  from  the  amplitude  samples  of  the  4FR  system. 

The  central  moments  pi .  pj,  .  .  .,  p„  of  a  random  variable  which  is  exponentially  or 
lognormal  distributed  are  well  known  and  can  be  used  for  comparison.  The  central  mo¬ 
ments  for  Uu?  logarithm  of  a  random  variable  which  is  exponentially  distributed  are  related 
to  Uu?  PoUgamma  function  (15).  The  Poligamma  function  'pb‘-l)(jr)  u  defined  as 

♦<“-»>(*)  -  |*„r(x  ♦  1J1  .  (31) 

dxv 


I 


where  F(x  ♦  1 )  is  Uie  Gamma  function.  The  central  moments  of  Uie  natural  iogariUun  of 
a  random  variable  which  is  exponentially  distributed  are  constants  given  by 
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from  the  exponential  distribution  (with  Rayleigh  envelope) 
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Fig.  13— Maximum  deviation  of  the  cumulative  distribution  of  terrain 
clutter  from  the  lognormal  distribution 
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and  'f'(Q)  *  0. §7721 57  is  the  difference  between  the  natural  logarithm  of  the  mean  value 
of  the  exponentially  distributed  random  variable  and  the  mean  value  of  the  natural  lo^hlhm 
of  the  random  variable.  The  numerical  values  of  the  !*oligumma  function  were  obtained 
from  Ibvts  (16). 

Similarly,  the  central  moment*  of  a  random  variable  (exponentially  distributed  in  am- 
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where  c  -  10  logio  e  with  e  -  2.7182818,  and  the  difference  between  the  median  and  the 
mean  value  of  the  random  variable  in  decibels  is  approximately  0.91. 

The  central  moments  of  ie  logarithm  of  a  random  variable  which  normally  is  log¬ 
normal  distributed  are  those  of  the  Gaussian  distribution 

n/2 

Mn  =  &l2J~J(2k-l),  (34) 

fe=l 

where  n  >  2  and  is  even.  If  n  is  odd,  nn  =  0. 

In  Figs.  16-55,  the  result  of  computing  the  central  moments  of  terrain  clutter  for  the 
four  frequencies  and  both  polarizations  are  shown.  In  these  plots  the  corresponding 
values  for  the  central  moments  for  the  exponential  and  the  lognormal  distribution  are 
shown  for  comparison.  Obviously,  these  results  show  that  terrain  clutter  in  general  is 
neither  exponential  nor  lognormal  distributed,  thus  supporting  the  results  of  the  Kolmogorov- 
Smimov  test  on  the  cumulative  distribution. 

One  obvious  conclusion  is  that  the  type  of  terrain  is  principally  reflected  in  the  mean 
value  of  the  clutter.  In  comparison  with  the  previous  empirical  results  obtained  for  sea 
clutter  (17),  now  we  find  a  greater  spread  in  all  the  moments. 


RESULTS  AND  CONCLUSIONS 

In  the  first  part  of  this  investigation  the  point-sc  terer  formulation  for  the  statistics 
o *  terrain  clutter  was  developed.  Although  the  model  is  quite  general,  specific  results  are 
given  only  for  the  case  of  a  side-looking  radar.  We  have  shown  that  many  of  the  statistical 
properties  of  terrain  clutter  can  be  derived  with  this  formulation.  An  important  result  of 
this  investigation  is  the  relatively  slow  rate  of  convergence  of  the  statistics  for  the  sum  of 
several  scintillating  signals  toward  the  Gaussian  distribution  as  the  number  of  components 
increases.  This  happens  because  the  mean  value  of  a  scintillating  signal  in  general  is 
nonzero. 


The  limiting  distributions  of  terrain  clutter  are  shown  to  be  the  Gaussian  distribution, 
that  obtained  over  homogeneous  terrains  such  as  deserts  and  farmland  for  example,  and  the 
distribution  of  a  scintillating  signal  produced  by  a  point-scatterer.  Terrain  clutter  which  has 
a  distribution  with  long  tails,  typical  of  lognormal  distributions,  should  be  encountered 
over  urban  areas  and  over  mountainous  terrain.  Figure  8  shows  selected  distributions  of 
terrain  clutter  obtained  with  the  4FR  system  which  have  interesting  features  tending  to 
support  the  point-scatterer  formulation  for  the  statistics  of  terrain  clutter. 


The  empirical  investigation  of  the  statistics  of  the  NRC  of  terrain,  in  decibel  units, 
in  general  has  neither  exponential  (of  Rayleigh  envelope)  nor  lognormal  distribution.  How¬ 
ever,  terrain  clutter  data  from  homogeneous  terrains  are  very  likely  to  be  exponentially 
distributed  (with  Rayleigh  envelope). 


Strictly  speaking,  terrain  clutter  is  never  lognormal  distributed  between  the  0  and  100 
percentiles.  However,  in  many  instances  the  distribution  of  terrain  has  long  tails  typical 
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of  lognormal  distributions,  and  for  these  cases  the  tail  of  the  distributions  may  be  approxi¬ 
mated  by  a  lognormal  distribution. 

Of  course,  in  order  to  evaluate  fully  the  merits  of  the  point-scatterer  formulation  for 
the  statistics  of  terrain  clutter  a  great  deal  of  further  investigation  of  the  statistics  of  data 
is  required. 
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